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INTRODUCTION 

The  Hyperbolic  Programming  Problem 

In  general,  the  role  of  mathematical  programming  can  be  summed  up  as  follows:  We 
are  looking  for  a  nonnegative  solution  of  the  variables  u^  U2,  . . . ,  u^  which  satisfies  the 
inequalities 


‘ra(“r-“2---"n)  ^ 

and  which  maximizes  the  given  function  v  =  v^Uj,  U2  . . . 

A  special  case  of  this  problem  is  known  as  linear  programming,  where  the  functions 
ff  . . .  f^  and  V  are  linear.  In  this  situation  the  given  inequalities  are 
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which  can  be  written  as  Uja^  +  U2a2  +  .  .  .  u^a^  <  a^,  or  In  the  following  vector  inequality 

foi  m  Au  s  a^.  The  function  v  can  be  written  as  v  =  CjUj  +  .  .  .  +  c^^u^  -  c^  or  v  =  c*u  -  c^. 

For  the  solution  of  this  problem,  as  it  is  known,  D^tzlg  discovered  the  simplex  method, 

which  we  presume  to  be  known  in  what  follows  [2,3]. 

This  study  discusses  another  special  case  of  mathematical  programming,  namely,  when 

the  functions  f*  ■  •  .  f„  are  linear  and  tne  function  v  is  a  linear  fractional  function.  Since  the 
1  m 

graph  of 


cx  -  c„ 

y  = - 

dx-do 

which  is  a  one-variable  linear  fractional  function,  is  hyperbolic  on  the  (x,  y)  coordinate  sys¬ 
tem,  we  call  this  a  hyperbolic  programming  problem. 

We  will  show  that  the  problem  of  hyperbolic  programming  can  be  solved  with  the  aid  of 
a  slightly  modified  simplex  method. 

The  method,  to  be  explained,  seems  useful  in  the  solution  of  economic  problems  where 
the  different  economical  activities  utilize  fixed  resources  in  propoi  tion  to  the  level  of  their 
values;  the  purpose  of  optimization,  however,  is  not  the  definition  of  a  revenue,  allocation  or 
economizing  extremum  (as  it  is  with  linear  programming),  but  the  extremum  of  a  specific 
index  number,  and  usually  the  most  favorable  ratio  of  revenues  and  allocations. 

On  the  basis  of  this  the  problem  of  hyperbolic  programming  is  the  following:  We  are 
looking  for  the  vector 


u  = 


u„ 


L“nj 


which  maximizes  the  function 


(1) 


x(u)  _  c*u  -  Cq 
y(u)  d*u  -  d^ 


subject  to  the  constraints  Au  ^  a^  and  u  li:  0.  The  set  of  those  points  that  satisfies  the 
constraints  we  will  call  the  set  of  feasible  solutions  and  designate  it  as  L.  It  is  well 
known  that  L  is  a  closed  convex  set  with  a  finite  number  of  extreme  points. 

Preliminary  Observations 

1.  The  necessary  and  sufficient  conditions  to  reduce  a  hyperbolic  programming  prob¬ 
lem  to  a  linear  programming  problem  could  be  easily  defined.  (I.e.,  the  points  of  set  R,  to  be 
defined  later,  should  fall  on  the  same  line.)  In  practice,  however,  we  do  not  think  it  necessary 
to  exclude  beforehand  all  degenerating  cases,  but  we  will  exclude  the  following  two  trivial 
cases,  (2)  and  (3): 
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The  set  L  such  that 


(2)  y  =  d*u  •  =  const. 

This  is  a  linear  programming  problem.  (If  y  ■  0,  then  the  problem  has  no  meaning.) 
The  set  L  such  that 


(3)  x(u)  =  v^y(u), 

where  =  const.  Thus  on  the  entire  set  of  L,  v(u)  =  except  for  the  common  0  points  of 
x(u)  and  y(u)  where,  however,  the  function  Is  not  Interpreted.  In  this  case  the  programming 
problem  has  no  meaning. 

2.  When  the  problem  Is  not  to  find  the  maximum,  but  the  minimum  of  the  function  v, 
we  may  equivalently  consider  the  problem  of  maximization  of  (-v). 

3.  We  will  assume,  throughout,  that  there  exists  at  least  one  point  In  the  set  L  such 
that  y(u)  is  positive.  This  can  be  done  without  loss  of  generality  since  we  may  always  con¬ 
sider  the  function  v  =  -x(u)/-y(u)  Instead  of  v  ::  x(u)/y(u).  (We  have  already  excluded  the  case 
y  8  0  In  (2).) 

Geometrical  Interpretation 

The  problem  can  also  be  formulated  as  follows: 

Let  us  consider  the  following  transformation; 


The  transformation,  as  can  be  seen,  consists  of  a  linear  transformation  and  a  parallel  shift. 
The  transformation  In  (4)  maps 


L  =  {u!Au  s  a^,  u  2  0} 

a  closed,  convex  set  with  finite  extreme  points  to 

R  =  {w(u)>  u(L } 

along  the  set  of  the  (y,  x)  plane.  (As  an  exception  we  choose  the  y  axis  to  be  horizontal.)  We 
can  say  the  following  about  the  set  R: 

1.  It  is  closed  and  convex; 

2.  Each  extreme  point  of  the  set  R  is  an  image  of  at  least  one  extreme  point  of  the 

set  L; 

3.  It  has  a  finite  number  of  extreme  points; 

4.  If  L  is  convex  polyhedron  then  R  Is  a  convex  polygon;  and 

0.  It  has  a  point  on  the  open  half  plane  y  >  0  (by  "Preliminary  Observations"). 
Statements  1  and  2  are  derived  directly  from  the  known  properties  of  linear  transformations, 
and  3  and  4  from  the  previous  statementr. 
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THE  SIMPLE  CASE 


Let  us  consider  the  (y,x)  plane.  The  value 
of  the  function  v(Uq)  =  x(’>Q)/y(UQ)  for  some  point 
u^  in  the  set  L  can  be  found  as  the  tangent  of  the 
angle  which  is  constructed  by  drawing  a  half 
radius  from  the  origin  to  the  point  w^ ,  which  in 
turn  is  the  image  point  of  some  point  u^ .  (See 
Figure  1.  The  dotted  line  indicates  the  case  when 
R  is  bounded.)  Thus  the  problem  is  to  find  such  a 
point  of  the  set  R ,  the  direction  tangent  of  which  is 
maximal,  or  to  find  such  points  of  set  L  which  were 
transformed  to  such  m:  xlmal  points  of  set  R.* 


The  Conditions  of  the  Simple  Case 

To  simplify  further  discussion,  we  will  present  a  simple  case  where  the  essence  of  the 
method  to  be  applied  may  be  observed. ' 

We  will  call  the  hyperbolic  programming  problem  a  simple  case  if  it  satisfies  the  fol¬ 
lowing  two  conditions; 

1.  The  set  of  possible  solutions  is  bounded  (therefore  it  is  a  convex  polyhedral  set, 
designated  as  P);  and 

2.  The  denominator  will  not  be  0  on  the  set  P.  This  constraint  along  with  the  assump¬ 
tion  in  the  second  section  (Preliminary  Observations)  —  considering  the  continuous  function  y  — 
means  that  on  the  entire  set  P , 


(5) 


y  =  d*u  -  d^  •  0. 


In  the  case  of  economical  programming  problem,  the  first  condition  means  thai  none  of  the 
economical  activities  to  be  programmed  can  be  unbounded.  The  assumption  in  the  second  con¬ 
dition  excludes  the  possibility  of  the  optimal  value  of  the  program  becoming  infinite.  It  can 
be  seen  that  these  conditions  correspond  to  practical  assumptions.  Even  for  this  reason  alone 
the  discussion  of  the  simple  problem  has  a  practical  Interest. 

Theorem  Concerning  the  Maximum  of  Function  v  for  the  Simple  Case 

With  the  aid  of  the  simplex  method,  we  can  examine  the  value  of  a  function  on  the  verti¬ 
ces  of  a  convex  polyhedral  set.  By  the  following  theorem  it  is  possible  to  apply  the  simplex 
method  to  the  simple  hyperbolic  programming  problem. 


•There  ih  .mother  poshible  ^eoniet  ric  .t  1  inte  rpret  .itlon .  Let  o  h  observe  the  f.iimly  of  hyper¬ 
planes  which  IS  delineii  by 

(c,  -  •  H,)u,  ♦  (c^  -  <1^)0^  ♦  .  .  .  ♦  -  'd^^, 

for  different  values  of  the  parameter  We  .ire  looking  lor  an  element  ol  the  f.imily  of  hyper¬ 
planes,  .It  which  the  value  of  the  p.irameter  ■  is  both  maximal  .ind  has  a  common  point  with 
the  set  L.  Obviously,  the  optim.il  hyperplane  cont.iins  an  extreme  point  from  the  set  L.  The 
difference  between  the  usual  geometrical  interpretation  of  line.ir  programming  and  this  c.ise 
is  that  in  the  former  trie  family  of  hyperplanes  is  parallel  and  in  the  latter  it  is  not.  This  in¬ 
terpretation  of  the  problem  was  worked  out  by  ,*.ndr,is  Prekop.i.  We  will,  however,  employ 
the  two  dimensional  transformation  for  the  sake  of  simplicity  in  what  follows. 

'We  consider  it  probable  that  the  majority  of  the  practical  -conomical  problems  will  satisly 
those  conditions  which  will  be  applied  in  this  case. 
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THEOREM  1:  If  the  linear  fractional  functional 


v(u) 


C*Uj^ 

d*u  - 

o 


has  the  following  properties: 

1.  the  area  of  definition  is  a  convex  polyhedral  P  defined  by  an  Au  ^  a^,  u  ^  0  con 
straint;  and 

2.  on  the  set  P , 


y  =  d*u  -  d^  >  0 , 

then  the  function  v  has  a  finite  maximum  on  the  set  P  which  is  achieved  on  at  least  one  vertex 
of  the  polyhedron. 

PROOF  OF  THEOREM  1:  Let  Pj,  P2  •  •  •  Pg  be  vertices  of  the  polyhedral  P.  Then, 
as  it  is  known,  any  z  of  P  can  be  written  as  a  convex  linear  combination  of  these  vectors, 

(6)  z  =*^iPi  +  ^2P2^**--^^V8 
where 

+  +  =  1 
and 

^  0,  i  =  1,  2,  ...  s. 

Let  us  assume  that  amcxig  the  vertices,  the  value  of  the  function  v  is  the  greatest  at 
Pfa,  i.e., 

(7)  v(pj,)  -  v(p,)  i  =  1,  2,  . . .  s, 
that  is 

x(Ph)  y(Pt)  ■  y(Ph)  • 

Multiply  both  sides  by  and  add  from  i  =  1  to  1  =  s 

(8)  x(pij)  l't»lty(Pi)  -  y(Ph)  *  (Pj)  . 


According  to  the  following  representation  and  by  considering  (6), 


d*z  -  dg  =  y(z) ; 


similarly. 


X  (pj)  =  x(z) . 
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Therefore,  from  (8) 


x(ph)  y(2)  -  y(Ph)  *  (*)i 


from  which  follows 

v^pjj)  ^  v(z) . 

With  this  we  proved  our  theorem  and  also  the  fact  that  with  the  aid  of  the  simplex 
method  there  exists  the  possibility  of  solving  a  simple  hyperbolic  programming  problem. 

The  Solution  of  a  Simple  Problem  by  the  Simplex  Method 

We  want  to  have  an  algorithm  with  the  aid  of  which  the  number  of  necessary  steps  In 
the  simplex  method  can  be  decreased  and  the  cyclic  Iteration  can  be  avoided. 

In  what  follows,  we  assume  that  the  Inequalities  have  already  been  converted  to  equations 
by  the  known  method  of  Introducing  slack  variables.  Then  we  can  write  the  equations  In  this 
form  Au  =  a^.  Let  uo  assume  that  we  know  a  basic  solution  of  the  problem  containing  m  positive 
components.  Let  us  further  assume  that  exactly  the  first  m  components  are  positive.  Then  the 
basic  vectors  are  a^,  aj,  ....  a^. 

^or  each  a^  (1  >  m  >  1,  m  -f  2,  . . . ,  n)  let  the  coordinates  In  terms  of  the  basic  vectors 
be  bjj ,  b2j ,  . . . ,  b^^^j  and  for  a^  the  coordinates  are  »  “2 »  *  ‘  »  “m  ’ 

Let 

*1  *  ‘^1^11  ♦  •••  ♦  ‘=m*^ml  ’  ‘^l' 

yj  =  djbn+...+d^b^j-dj, 

(1  =  m  4  1,  m  4^  2,  . . . ,  n),  and 
=  xyi  -  yxj . 

(If  aj  Is  a  basis-vector,  then  bjj  =  0»  j ;  bjj  =  1,  and  accordingly  Xj  =  y^  =  tj  =  0.) 

Now  we  can  write  the  following  simplex  tableau  which  is  an  extended  one  compared  to 
the  one  used  In  linear  programming; 


<^1- 

"1 

C2  .  .  . 

m 

•^m+l  •  •  • 

‘^k 

.  .  .  c„ 

^o 

‘‘i 

d2  .  .  . 

‘^m 

‘^m+l  •  •  • 

‘Ik 

•  •  •  ‘*n 

basis 

»0 

*1 

®2  •  •  • 

»m 

*m+l 

.  .  .  a 

n 

"1 

‘•l 

*1 

“1 

1 

0  ... 

0 

^l,m4l  •  •  • 

‘’Ik 

•  •  •  ‘’in 

^2 

<*2 

*2 

“2 

0 

1  ... 

0 

bo  •  •  • 

2,m+ 

‘’2k 

•  •  •  ‘’2n 

, 

‘‘m 

*m 

“m 

0 

0  ... 

1 

^m,m+l  ■ 

‘’mk 

‘’mn 

X 

0 

0  ... 

0 

*m4l  •  •  • 

’Sc 

.  .  .  Xn 

y 

0 

0  ... 

0 

^m+1 

^k 

•  •  •  ^n 

‘r 

x/y 

0 

0  ... 

0 

%+l  •  • 

‘k 

•  •  •  ‘n 
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The  arrangement  of  the  known  data  is  understood  from  the  tableau.  Into  row  t^  of  col- 
unm  a^  we  write  the  function  value  v  =  x/y  given  by  the  current  basic  solution. 

Let  us  view  the  positive  coordinates  of  the  vector  a^^  and  let 

min  —  =  6  . 

We  have  assumed  that  Uj  >  0,  therefore  6  >  0.  (We  will  return  to  the  question  of  degeneracy.) 
As  is  known,  the  value  of  the  linear  functions  x(u)  and  y(u)  becomes 

X'  =  X  - 

y'  =  y  - 


after  the  vector  a^^  has  been  brought  into  the  basis. 

Thus  by  exchanging  the  basis-vector  the  changed  value  of  the  function  v  is: 

,  '■'’Sc 

v'  =  -  . 

y-% 

v'  becomes  greater  than  v,  l.e.,  the  changing  to  the  new  basis  means  approaching  to  the  maxi¬ 
mum  if 


or 


y  - 


0. 


y(y  -  sy,,) 


>  0. 


Considering  that  in  the  case  examined  y  -  0,  y  -  dy^^  ^  0,  and  6  ^  0,  we  can  say  the  following: 
If  we  bring  into  the  basis  a  vector  a^  having  a  positive  coordinate  for  which  the  t^  la  positive 
in  the  last  row  of  the  simplex  tableau,  t^  =  xy^  ~  yx^  ^  0,  then  the  value  of  the  function  v  will 
Increase. 

Therefore  the  simplex  method  employed  in  hyperbolic  programming  consists  of  bring¬ 
ing  into  the  basis  successively  all  positive  t^  vectors  (if  they  have  positive  coordinates).  The 
method  of  calculating  all  u^ ,  b^^ ,  x  ,  y  ,  x^ ,  y^  values  corresponds  to  the  one  used  in  linear 
programming;  only  v  =  x/y  and  t^  =  xy^  -  yx^  have  to  be  computed  by  using  these  special 
formulas. 
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Since  (1)  the  number  of  vertex  points  are  finite,  (2)  the  bringing  in  of  each  vector  char¬ 
acterized  by  t^  >  0  increases  the  value  of  v,  and  (3)  the  function  v  will  reach  a  maximum 
value  at  a  certain  vertex  point,  we  must  achieve  an  optimal  solution  after  a  finite  number  of 
steps  by  bringing  the  positive  t^  vectors  into  the  basis. 

The  tableau  containing  an  optimal  solution  is  an  optimal  tableau.  From  the  above  it 
follows  that  in  the  optinul  tableau  all  t^  =  0.  During  the  discussion  of  the  general  case  (section 
entitled  "Application  of  the  Simplex  Method  to  the  Solution  of  the  'Good'  Case  of  the  General 
Problem,"  tableau  2)  we  will  prove  that  this  condition  is  sufficient:  Each  such  tableau  that 
contains  t^  -  0  is  an  optimal  tableau. 

The  number  of  computations  necessary  to  solve  a  simple  hyperbolic  programming  prob¬ 
lem  with  the  aid  of  the  simplex  method  exceeds  only  by  a  small  percentage  the  number  of  oper¬ 
ations  required  for  a  linear  programming  problem  of  the  same  size.  Programs  for  digital 
computers  can  be  applied,  with  a  slight  alteration,  to  hyperbolic  programming. 

THE  GENERAL  PROBLEM 
Definitions  —  Existence  Theorems 

In  the  case  of  the  general  problem  we  do  not  assume  that  the  set  L  is  bounded  or  that 
the  denominator  has  no  0  point  on  the  set  L.  What  makes  the  discussion  of  the  general  problem 
necessary,  as  we  will  later  see,  is  the  fact  that  in  certain  cases  where  the  conditions  for  the 
simple  problem  are  absent  the  problem  may  still  have  an  optimal  solution  even  though  the  set 
L  is  unbounded  or  the  denominator  becomes  0. 

Definitions. 

1.  We  will  call  the  "good"  point  of  the  set  L  where  either  y{u)  >  0  or  y(u)  =  0,  x(u)  <  0. 
(According  to  the  third  observation  under  "Preliminary  Observations,"  L  always  has  a  "good" 
point.) 

2.  We  will  call  the  "bad"  point  of  the  set  L  where  either  y(u)  <  0  or  y(u)  »  0,  x(u)  ^  0. 

3.  The  singular  point  of  the  set  L  is  at  y(u)  =  x(u)  =  0. 

Therefore,  in  the  transformation  of  the  third  section  (Geometrical  Interpretation),  the  good 
points  are  mapped  on  the  inside  of  quadrants  I  and  IV,  on  the  positive  side  of  the  y  axis  (hori¬ 
zontal),  and  on  the  negative  side  of  the  x  axis;  the  bad  points  are  mapped  on  the  inside  of 
regions  n  undin,  on  the  negative  part  of  the  y,  and  the  positive  part  of  the  x  axis;  the  singular 
points  are  mapped  into  the  origin. 

We  will  prove  the  following  theorems: 

THEOREM  2  ("BAD"  CASE):  If  the  set  L  has  bad  points  then  the  function  v  is  not 
bounded  from  above,  and  the  programming  problem  has  no  optimal  solution.  (See  Figures  2 
and  3.) 


THEOREM  3  ("GOOD"  CASE):  If  the  set  L  has  only  good  points  and  the  function  v 
has  a  finite  maximum  on  the  set  L,  then  this  is  taken  on  at  least  one  extreme  point  of  the  set 
L.  (See  Figures  4  and  5.) 

THEOREM  4  (SINGULAR  CASE):  If  the  set  L  has  no  bad  points,  but  has  singular 
points  (and  naturally  good  points)  then  it  has  a  good  point  where  the  function  v  takes  on  a  finite 
maximum  value.  (See  Figures  6  and  7.) 
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Figure  5  Figure  6  Figure  7 


Proof  of  the  Existence  Theorems 

PROOF  OF  THEOREM  2  (BAD  CASE):  On  the  basis  of  the  statement  supporting  Eq. 
(3),  for  each  bad  point  (P^.)  of  the  set  L  we  can  choose  such  a  good  point  (Pj)  that  on  the  line 
segment  Pj.Pj  connecting  these  two  points  there  will  be  no  singular  points.  Therefore  the  line 

segment  PJ.Pj  (the  image  on  the  plane  of  Pj)  has  a  P’  common  point  with  the  x  axis  which 
is  different  from  the  origin.  If  P'  is  on  the  positive  x  axis,  then  in  the  quadrant  I,  if  it  is  on 
the  negative  x  axis  then  in  the  quadrant  III  we  can  choose  along  the  line  segment  P|.Pj  such  a 
series  of  pomts  converging  to  P',  along  which  series  the  direction  tangent  grows  through  posi¬ 
tive  numbers  beyond  ail  limits.  Therefore  tlie  f'jnction  v  is  aot  bounded  from  above,  and  the 
programming  problem  has  no  solution. 

The  theorem,  of  course,  can  be  proved  without  the  geometrical  interpretation;  however, 
we  will  exclude  this  method  of  proof. 

PRCX)F  OF  THEOREM  3  (GOOD  CASE):  The  theorem  can  be  considered  as  the  gen¬ 
eralization  of  Theorem  1.  Since  if  we  do  not  declare  the  set  L  bounded,  then  the  function  v 
does  not  always  have  a  maximum  on  the  set  L.  Furthermore,  as  we  will  see,  the  function  v, 
bounded  from  above,  is  not  a  sufficient  condition  for  the  existence  of  the  maximum. 

Those  points  of  the  set  L  for  which  y  =  0,  x  <.  0  cannot  be  optimal.  These  points  we 
can  approach  only  by  such  points  the  images  of  which  are  in  quadrant  IV.  This  way,  however, 
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the  function  v  becomes  infinite  only  through  negative  values,  therefore  It  does  not  reach  a 
maximum  value.  Consequently,  it  Is  sufficient  to  prove  Theorem  3  for  those  points  only  for 
which 

(9)  y(u)  »  d*u  -  d^  >  0. 

PROOF:  Let  Pp  P2,  •  •  ■  i  Pg  1>^  the  extreme  points  of  the  set  L  and  the  value  of  v  be 
the  highest  at  Pj^,  l.e., 

(10)  ''(Ph)  - ’'(Pl)*  1*1,2,  ...s. 

We  prove  that  if  the  set  L  has  a  point  z  where  th  i  function  v  is  greater  than  at  the 
chosen  Pj^  extreme  point,  l.e., 

(11)  v(z)  >v(pj^), 
then  it  also  has  such  a  point  z  where 

(12)  v(z)>v(z). 

Obviously  z  cannot  be  an  extreme  point  of  the  set  L.  It  follows  that  we  can  choose  a  Pg 

extreme  point  to  point  z,  and  a  number  o  >  0  so  that  if  q  =  z  •  p^,  then  z  =  z  -f  oq 

=  (1  ^a)z  -  op  should  still  be  a  point  of  the  set  L;  however, 

(13)  v(z)  =  -  -  -  =  -  . 

d-i-d„  'l.<.)d*2-»d*Pg-d„  (l.a)y(2)-o,(Pg) 

From  (11)  and  (10)  it  follows  that  v(z)  >  v^Pg),  i.e.,  x(z)y^Pgj  >  x^p^jyfz).  When  both  sides 
are  multiplied  by  {-a)  and  when  (1 -f of)x(z) y(z)  is  added  we  obtain 

(1 +a)x(z)  -  ox^Pgj  . 

Equation  (12)  follows  from  this  by  considering  (13).  Thus  the  theorem  is  proved. 

As  far  as  the  maximum  of  the  function  v  is  concerned  the  theorem  permits  three  cases: 

1.  The  function  v  reaches  its  maximum  on  the  set  L  and  it  is  reached  on  some  ex¬ 
treme  point  of  the  set  L  (Figure  4). 

2.  The  fimctlon  v  is  not  bounded  on  the  set  L  from  above  (Figure  8). 

3.  The  function  v  is  bounded  from  above  on  the  set  L,  but  it  does  not  reach  its  maxi¬ 

mum  (Figure  9). 

These  cases  will  be  differentiated  by  examining  the  extreme  points  in  a  way  to  be  given 
in  the  section  entitled  "Application  of  the  Simplex  Method  to  the  Solution  of  the  'Good'  Case  of 
the  General  Problem." 


x(z) 


(l  +  o)y(z)  -  oy  ^Pgj  <  y(z) 
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PROOF  OF  THEOREM  4  (SINGULAR  CASE): 

This  theorem  will  be  proved  by  the  presentation  of  the  optimal  solution  in  the  section 
entitled  "Solution  of  the  Singular  Case." 

Preparation  of  the  Solution:  Differentiation  of  the  Cases 

The  hyperbolic  programming  problem  solved  by  the  simplex  method  can  be  initiated  in 
the  following  steps:  * 

1.  We  define  a  starting  solution  by  the  same  method  that  is  used  l^  linear  programming. 
(The  inequalities  are  converted  to  equations  by  the  introduction  of  slack  variables,  and  so  on.) 

2.  We  examine  the  y  value  belonging  to  the  starting  solution  —  (a)  if  y  is  positive  we 
proceed  to  step  3;  (b)  if  y  is  negative  we  proceed  to  functions  (-x)  and  (-y)  instead  of  func¬ 
tions  X  and  y  and  continue  with  step  3;  and  (c)  if  y  is  0  we  see  if  the  value  of  y  could  be  in¬ 
creased  by  bringing  any  vector  into  the  basis.  If  it  can  be  Increased  we  proceed  according  to 
(a);  if  it  cannot  be  Increased,  only  decreased,  we  act  according  to  fb);  if  it  cannot  be  either  in¬ 
creased  or  decreased  it  corresponds  to  the  excluded  case  y  i  0. 

3.  With  linear  programming  we  start  searching  for  the  minimum  of  the  function  y  — 

(a)  if  in  the  meanwhile  we  reach  a  negative  y  value  we  Interrupt  the  procedure,  there  is  no 
solution  (Bad  case);  (b)  if  min  y  ■  0  then  proceed  to  step  5  (Good  case);  and  (c)  if  min  y  =  0 
then  proceed  to  step  4. 

4.  We  add  to  the  original  given  equations  another  condition. 


and  find,  with  linear  programming,  the  maximum  of  the  function  x  taking  the  added  conditio.. s 
into  consideration'  —  (a)  if,  in  the  meanwhile,  we  reach  a  positive  x  value  we  Interrupt  the 
procedure,  there  is  no  solution  (Bad  case);  (b)  if  max  x  <  0,  then  proceed  to  step  5  (Good  case); 
and  (c)  if  max  x  =  0,  then  proceed  to  the  procedure  to  be  described  in  the  section  entitled  "So¬ 
lution  of  the  Singular  Case." 


*11  know  I  roll)  the  c  li.i  r.tctc  r  of  the  problem  th.tt  min  y  ■  0  then  these  initi.itin^  ste[)s  .ire  not 
needed 

I.e.,  lA  e  .ire  brin^in^  those  .t|^  vectors  into  the  basiN  for  which  yj^  U,  U. 
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5.  In  the  course  of  the  linear  programming  procedure,  performed  In  step  3,  we  com¬ 
pute  the  corresponding  v  =  x/y  value  for  each  tableau.  If  3(b)  or  4(b)  were  the  case,  It  would 
be  best  to  choose  as  an  Initial  tableau  the  one  corresponding  to  the  highest  v  value. 

In  case  of  a  more  complex  problem,  steps  3  and  4  should  be  performed  by  a  computer. 

Application  of  the  Simplex  Method  to  the  Solution  of  the 
"Good"  Case  of  the  General  Problem 

It  follows  from  Theorem  3  that  In  order  to  obtain  an  optimal  point  it  is  sufficient  to 
examine  the  extreme  points  of  the  set  L.  Therefore  the  simplex  method  can  be  applied. 

We  remind  the  reader  that  In  the  section  entitled  "The  Solution  of  a  Simple  Problem  by 
the  Simplex  Method,"  we  did  not  use  the  assumption  that  the  set  L  is  bounded  and  furthermore 
we  needed  only  the  milder  assumption  y  -  0  Instead  of  y  0. 

This  assumption  (y  =  0),  however,  holds  for  the  "good"  case  of  the  general  problem. 
Accordingly,  the  Simplex  Method  can  be  applied  without  alteration.  What  we  do  not  know  at  this 
point  is  whether  the  problem  has  an  optimal  solution.  Through  the  procedure,  however,  we  will 
find  the  answer. 

In  the  course  of  the  procedure  we  may  obtain  the  following  tableaux; 

T.I.  In  the  last  row  of  the  tableau  there  is  a  t^  >0  such  that  the  corresponding  aj^ 
vectors  have  positive  coordinates.  In  this  case  the  given  tableau  is  not  optimal  and  the  proce¬ 
dure  is  continued  by  bringing  vector  a^^  into  the  basis. 

T.2.  There  is  no  positive  number  m  the  last  row  of  the  tableau,  all  t.  •  0.  In  this  case 
the  function  v  reaches  its  maximum  on  the  set  L  and  the  solution  is  optimal. 

T.3.  There  is  a  t^  <  0  in  the  last  row  of  the  tableau,  but  one  of  the  corresponding  aj^ 
vectors  has  no  positive  coordinate.  In  this  case  the  function  v  does  not  reach  the  maximum  on 
the  set  L,  and  the  programming  problem  has  no  solution.  We  Interrupt  the  procedure. 

As  can  be  seen,  these  three  cases  are  identical  with  the  ones  in  linear  programming. 

The  only  difference  is  that  in  the  case  of  T.3  the  function  v  can  be  bounded  from  above,  while 
the  corresponding  linear  function  cannot  be.  (See  the  discussion  of  the  case  T.3,  below,  and 
Figure  9.) 

Proof  for  the  T.I  Case.  Our  statement  about  the  T.I  case  does  not  need  to  be  proved 
according  to  the  section  on  tue  simplex  method  of  solving  simple  problems.  From  the  finite 
extreme  points  of  the  set  L  it  follows  that  the  T.I  type  tableau  may  occur  only  in  a  finite  num¬ 
ber  of  steps;  and  afterwards  we  must  arrive  at  a  T.2  or  T.3  type  tableau. 

Proof  for  the  T.2  Case.  Let  us  assume  that 

“■  . “n.’" . ”]■ 

a  possible  solution  appears  in  the  tableau,  i.e., 

m 

(14)  “j  “J  =  “o  ■ 

)■! 

and  there  is  no  positive  number  in  the  last  row  of  the  tableau,  i.e.. 


(15) 


t.  =  x(u)yj  -  y(u)Xj  0  (i  =  1,2 . n). 
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We  prove  that  if  z*  =  Z2,  .  . . ,  z^  j  is  a  possible  solution,  i.e. 


(16) 


11  "l“l  =  “o- 
1=1 


then 


(17) 


v(u)  -  v(z). 


Any  a^  vector  on  the  a^,  . . . ,  a^  basis  can  be  expressed  as 


m 


j  -  ^  ^ji^j’  (i  -  If  2, . . . , n), 


j=l 


inserted  into  (16)  this  ((ives 


n  m 


m  /  n 


11  Vi  -llhll  -  Hill  =“o- 

i=l  i=l  j=l  j=l  \i=l 


If  we  combine  this  with  (14)  (and  consider  that  a^  can  be  expressed  only  uniquely  on  the 
Uj . a^^  basis)  we  get 


(18) 


II  =  “)• 


i=l 


By  multiplying  (15)  by  z.  and  adding  from  i  =  1  to  i  =  n,  we  get 


(19) 


•• 

x(u)  22  22  ’‘i^ 


i  =  l 


i=l 


Considering  the  definition  of  y^  and  (18),  we  may  transform 


n 

rn 


i-i 


in  the  following  manner: 


n  n/m  \  m/n 

(Cont.)  i  =  l  i  =  l  \j  =  l  1  j  =  l  \i  =  l 


(20) 


II  \  II 

H  ^  = 


i  =  l 
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(20) 


y(u)  -  y{z) , 


Similarly, 


n 

(21)  ’‘i^i  ^  ‘ 

1=1 


By  using  (20)  and  (21)  and  substituting  Into  (19),  we  get 

x(u)  [y(u)  -  y(z)J  -  y(u)  tx(u)  -  x(z)J, 


that  Is, 

(22)  x(u)y(z)  -y{u)x{z). 

We  consider  that: 

1.  y(u)  =  0  Is  Impossible.  Since,  If  y(u)  =  0  then  x(u)  ^  0  and  choosing  a  negative 

tj^  =  :;(u)yj^  -  y(u)Xj^  =  x(u)yjj  <  0  from  which  follows  yj^  ^  0.  This  way,  however,  by  bringing 
aj^  into  the  basis  we  would  reach  a  bad  point,  but  the  set  does  not  have  a  bad  point. 

2.  If  y(u)  ^  0,  y(z)  =  0,  then  v(z)  =  -  «  and  thus  (17)  Is  true. 

3.  If  y(u)  *  0,  and  y(z)  0,  then  (17)  follows  from  (22). 

With  this  our  statement  about  the  case  T.2  has  been  proved.  This  proof  also  proves 
the  penultimate  paragraph  of  the  section  on  the  simplex  method  of  solving  simple  problems. 

Proof  for  the  T.3  Case.  To  maintain  the  assumption  in  (14),  let  us  assume  that  vector 
a^^  (k  -  m  1)  Is  the  one  for  which 

123)  t^  =  x(u)  yj^  -  y(u/  \  ^  0 , 

and  aj^.  In  terms  of  the  given  basis,  has  no  positive  coordinates,  that  is, 

m 

(24)  a^  =2^  ' 

j  =  l 

where 

(25)  0  =  1.2,...,m). 

Through  multiplying  Eq.  (24)  by  a  number  A  >  0  and  by  adding  Eq.  (14)  to  it,  we  obtain 

m 

L  ’“o- 

j=l 

By  considering  (25),  u.  -  Ab^j  ^  0,  and  thus  the  following. 
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7.* 


2k' 


%  -  ^‘^mk' 


Is  a  possible  solution.  The  positive  coordinates  of  this  may  take  on  any  large  value;  therefore 
the  set  L  is  not  bounded  In  this  case. 

Let  us  consider  the  function  x  evaluated  at  the  point  z: 


m 


:(z)  =  ‘=j“j  -  *"0  - 

j=l  j=l 


m 

S  " 

i _ _ 

J=1 


‘"jV 


=  x(u)  -  Ax^ 


Similarly  y(z)  expressed,  provides 


(26) 


x(u)  -  Ax. 

V(z)  =  - ^  . 

y(u)  -  Ay^ 


The  value  of  v(z)  Increases  with  the  Increase  of  A,  since  considering  (23) 

^  _  x(u)yk  •  yM\  ^  ^ 

Since 


y(z)  =  y(u)  -  Ayj^  ^  0, 

Inequality  must  hold  at  any  a  ■  0,  It  follows  that 

0: 


1.  If  yjj  =  0,  then 


x(u)  -  Ax. 

v(z)  = - ^ 

y(u) 


and  according  to  (23)  |^-y(u)Xj^j  •  0,  therefore  Xj^  <  0.  Thus  with  the  increase  of  A,  v(z)  In¬ 
creases  without  bound. 

2.  Uy,  ^  0,  then,  as  can  be  easily  shown. 


v(z)  <  —  =  lim  v(z) . 

A  •  +  X 

Therefore  the  function  v  is  bounded  from  above  on  the  set  L,  Its  lowest  bound  is  but 

It  does  not  reach  this  value  on  any  finite  point  of  the  set  L.  With  this  our  statements  for  the 
T.3  case  are  proved. 


♦If  aj^  18  the  only  vector  satisfying  the  assumptions  in  T.i.  If  there  are  more  such  agj,  •'k' 

vectors  then  the  least  upper  bound  is  max  ^kf^^kf  . ^ 


150 


BELA  MARTOS 


Solution  of  the  Singular  Case 

In  the  singular  case  the  set  L  has  singular  points  where  the  fractional  function  v  Is 
not  Interpreted,  since  both  Its  numerator  and  denominator  become  0.  From  a  practical  point 
of  view  this  kind  of  problem  should  be  considered  solvable  when  there  Is  such  a  point  among 
the  nonsingular  points  of  the  set  L  where  the  function  v  reaches  its  maximum  value.  We  stated 
that  If  the  set  L  has  no  bad  points  then  there  is  always  a  solution  In  this  sense.  (In  other  words, 
even  though  the  function  has  no  boundary  value  In  the  singular  point,  it  has  llm.  sup.) 

To  solve  the  singular  case  we  start  out  from  the  tableau  that  was  obtained  In  step  4(c) 
of  "Preparation  of  the  Solution:  Differentiation  of  the  Cases,"  In  which  y  =  0  and  x  =  0.  (Sin¬ 
gular  *''bleau.)  From  this  It  followo  that  In  the  last  row  of  the  singular  tableau  all  t^  =  0  can 
be  found.  We  have  to  return  to  the  y^  and  x^  values  when  examining  the  tableau: 

1.  The  singular  tableau,  considering  the  minimum  of  the  function  y,  is  an  optimal  one, 
therefore  each  y^  0.  (But  among  these  there  is  yj^  ^  0,  otherwise  there  were  y  =  0.) 

2.  The  singular  tableau  Is  an  optimal  one  for  the  maximum  problem  of  the  function  x 
which  Is  enlarged  bv  the  assumption  d*u  =  d^,  i.e.,  for  all  such  vectors  to  which  yj^  =  0  value 
belongs,  Xj^  -  0. 

There  are  two  kinds  of  singular  tableaux: 

5.1.  Some  vector,  having  ^  0  and  maximum  yj^  value  has  positive  coordinates. 
Bringing  this  vector  into  the  base  we  obtain  an  optimal  solution  immediately. 

5.2.  No  vector,  having  y^^  •  0  and  maximum  Xj^/y|^  value,  has  positive  coordinates. 

As  described  below,  we  can  again  obtain  an  optimal  (not  basic)  solution. 

Proof  for  the  S.l  Case.  We  proved  in  the  T.2  case  that  a  tableau  is  nonsingular  and 
optimal  If  (1)  y  -0  and  (2)  tj  •  0  (i  =  1,  2, . . . ,  n). 

1.  After  bringing  in  the  aj^  vector  acceptable  to  the  conditions  in  the  S.l  case,  the 
changed  y'  value  of  the  new  tableau  will  be  positive 

y'  =  y  -  <5y^  =  -  5y,^  '0, 

since  yj^  <  0.  Therefore  the  new  tableau  is  not  singular. 

2.  We  obtain  the  altered  tj  values  as  follows: 

f 

<  =  -  ,  where  g  is  the  index  of  the  vector  that  was  taken  out  of  the  basis. 

I  V 

t; .  x'y;  -  y'x; '  -  «Xk(y,  -  «,yk)  •  »yk(x,  -  ’>i\)  =  -  6(\yi  -  y|,X|). 

In  order  to  prove  that  t|  •  0  for  all  i,  we  have  to  show  that 

(27)  \yi-yk^=o 

for  all  i,  unless 

\  \ 

yk  ~  yi  ’ 


(28) 
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for  all  such  i  when  the  latter  fraction  Is  finite  and  not  undetermined  (i.e.,  if  we  brought  in  the 
vector  having  a  maximal  x^/yj). 

Regarding  the  statement  in  (27), 

1.  is  trivial  in  case  Yj  =  =  0| 

2.  case  =  0,  x^  0  was  excluded, 

3.  is  true  in  case  y^  =  0,  Xj  0  because  yj^  0,  and 

4.  in  case  y^  •  0,  x.  y^  is  finite  and  thus  (27)  follows  from  (28). 

With  this  we  proved  that  we  reached  an  optimal  tableau  in  the  given  manner. 

Proof  tor  the  S.2  Case.  Let  us  choose  a  vector  a^,  (y^^  <  0;  b^j^  0,  j  =  1,  2,  . . . ,  m) 
for  which  value  is  maximal  compared  to  all  finite  Xj/yj  values;  i.e.,  (28)  holds. 

Thus,  as  we  saw  in  the  proof  of  the  S.  1  case. 


(27)  \yi  '  yk^  i  =  1,2 . n, 

is  true  for  all  i.  Let  us  assume  that  the  basis  of  the  singular  tableau  happens  to  be  a^ 
a^^  and  the  solution  in  it  is 

s*  =  [Sj,  Sj,  ...,  8^,  0,  ...,  O]. 

Therefore, 


x(s)  =  y(8)  =  0. 


Let  us  now  formulate  the  solution 


(29) 


u*  = 


(k) 

-  ■^b||^,  Sn  -  Ab«|^,  . . . ,  8  -  Ab  1^,  0  . . . ,  A ,  . . .  0 


where  A  is  an  arbitrary  positive  number.  We  showed  in  the  proof  of  the  T.3  case  that  the 
solution  u  is  possible. 

On  the  example  of  (26)  we  can  see  that 


(30) 


X(s)  -  Ax.  X. 

v(u)  =  - ^  . 

y(8)  -  Ay^  y^ 


Now  we  will  show  that  the  solution  u  is  optimal  too.  If 

. ^n] 

is  an  arbitrary  possible  solution  for  which  y(z)  -  0,  then 
(31)  v(u)  -  v(z). 

We  multiply  both  sides  of  the  inequality  in  (27)  by  z.  and  add  from  i  =  I  to  i  =  n, 

n  n 

\  ^  Vi^-  yk  _  ’‘in- 
i=l  i=l 


(32) 
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n 


1=1 


z 


I 


can  be  changed,  on  the  example  of  (20),  to 

n 

=  y(8)  -  y(2)  =  -  y(z) . 

Similarly, 


n 

V" 

tn 


XiZj  = 


-  x(z)  . 


Thus  from  (32)  we  have 


-x^^y(z)  -  -  y,^x(z). 

Thus  taking  Into  consideration  that  yj^  0  and  y(z)  >  0 


\  _  x(z) 

y^  ~  y(z)  ’ 


Therefore,  considering  (30),  (31)  is  true  and  the  solution  of  (29)  is  optimal. 

We  thus  proved  our  statement  for  the  S.2  case  and  at  the  same  time  produced  the  method 
for  an  optimal  solution.  With  this  Theorem  4,  the  existence  theorem,  was  also  proved. 

Supplements 

In  this  section  we  wish  to  give  supplementary  background,  without  proof,  that  is  to  say 
we  wish  to  give  solutions  to  a  few  problems  that  arise  during  hyperbolic  programming. 

The  Case  of  De^ej^erac^.  The  degeneracy  depends  on  the  mutual  position  in  space  of 
the  polyhedron  generated  by  the  column  vectors  of  the  matrix  A  and  the  a^^  vector.  This  posi¬ 
tion  is  the  same  in  both  hyperbolic  and  linear  programming.  It  can  be  presumed  that  the 
degeneracy  problems  arising  in  the  course  of  hyperbolic  programming  can  usually  be  solved 
with  the  aid  of  the  perturbation  method  worked  out  by  Charnes  [ij.  We  have  to  call  attention, 
however,  to  the  fact  that  in  case  of  degeneracy  the  handling  of  singular  tableaux  becomes  more 
complicated  when  the  solution  of  the  problem  might  be  reached  only  by  computing  more  suc¬ 
cessive  singular  tableaux.  It  remains  to  be  examined  whether  in  this  case  the  perturbation 
method  is  appropriate  for  avoiding  cycling. 

The  Enlarged  Problem.  In  case  not  only 

a,,  u,  +  a.«  u„  4  . . .  +  a.  u„  -  a.  •  0 

il  1  12  2  m  n  lo 


^ri^^2“2^ 


in  n  lo 


inequalities,  but  also 
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equations  appear  in  the  original  conditions,  we  can  change  Au  ^  a^  into  Au  Et  =  a^.  This 
way  we  obtain  a  starting  solution  easily  and  ensure  the  linear  Independence  of  the  row  vectors 
of  the  enlarged  [AE]  matrix. 

This  way,  however,  the  problem  is  interpreted  on  an  enlarged  set  L'  rather  than  on  the 
original  set  L.  If,  originally,  there  were  equations  (and  not  inequalities)  in  rows  j  j2>  ■  •  •  j| 
then,  in  order  to  obtain  optimal  solution  on  the  original  L  only,  we  would  have  to  ensure  that 


be  in  the  optimal  solution.  For  this  reason  the  function  to  be  maximized  must  be  brought  into 
the  following  form: 


c,u,  .  C;..;  •  ■  ■  ■  •  -  Ntj,  -  .  .  .  -  Nt.,  - 

d,U,  .  djUj  .  .  .  .  .  d„u„  -  d^ 

where  N  is  a  very  large  positive  number.  As  in  linear  programming  there  is  no  need  for  the 
numerical  definition  of  N.  The  simplex  tableau,  however,  must  be  enlarged  by  two  rows  which 
contain  the  coefficient  of  N  In  the  expressions  x.  and  t^. 

All  Optimal  Solutions.  In  hyperbolic  programming  all  optimal  points  form  a  convex  set 
which  has  a  finite  number  of  extreme  points.  No  extreme  points  belong  to  the  set  of  optimal 
points  in  the  S.2  singular  case,  since  the  optimum  set  is  generated  partially  by  infinitely  dis¬ 
tant  points  and  partially  by  singular  points. 

The  method  for  the  definition  of  all  optimal  solutions  is  presented  only  for  the  case 
where  the  set  of  optimal  points  is  bounded.  If,  in  the  procedure  to  be  described,  we  reach  an 
optimal  tableau,  having  a  vector  for  which  t^  =  0  and  all  of  its  coordinates  are  not  positive, 
then  the  set  of  optimal  points  is  not  bounded.  This  criterion  holds  even  when  we  arrive  at  a 
singular  tableau  while  determining  all  optimal  points.  (In  the  S.2  case,  as  we  know,  the  optimal 
set  is  not  bounded.) 

If  the  above  described  case  does  not  present  itself  then  the  set  of  all  optimal  points  is 
bounded.  In  this  case  the  extreme  points  of  the  set  of  optimal  points  can  be  obtained  by  start¬ 
ing  from  an  optimal  tableau  and  successively  bringing  into  the  basis  all  those  vectors  for 
which  t^  =  0.  In  the  meanwhile  we  may  also  obtain  singular  tableaux  which  should  be  taken  into 
considerations  in  further  operations.  All  optimal  solutions  can  be  obtained  as  convex  linear 
combinations  of  solutions  from  all  tableaux  with  the  constraint  that  the  coefficients  of  solutions 
from  all  nonsingular  tableaux  cannot  be  zero. 

The  image  of  the  set  of  all  optimal  points  may  be: 

1.  an  extreme  point  of  the  set  R  (Figure  10); 

2.  the  closed  line  segment  connecting  two  neighboring  extreme  points  of  the  set  R 
(Figure  11); 

3.  a  half  radius  originating  from  one  extreme  point  of  the  set  R  (Figure  12); 

4.  a  half-open  line  segment  connecting  the  origin  with  an  extreme  point  of  the  set  R 
(Figure  6);  and 

5.  an  open  half  radius  starting  from  the  origin  (Figure  7). 

Cases  3  and  5  can  occur  only  when  the  set  L  is  not  bounded.  Cases  4  and  5  can  occur  if,  and 
only  if,  there  is  a  singular  one  among  the  tableaux. 
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Figure  10  Figure  II  Figure  1 
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